Does the Heisenberg model describe the muhimagnon spin dynamics in 
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We compute the absorption spectrum for multimagnon excitations assisted by phonons in insu- 
lating layered cuprates using exact diagonalization in clusters of up to 32 sites. The resulting line 
shape is very sensitive to the underlying magnetic Hamiltonian describing the spin dynamics. For 
the usual Heisenberg description of undoped Cu-O planes we find, in accordance with experiment, 
a two-magnon peak followed by high energy side bands. However the relative weight of the side 
bands is too small to reproduce the experiment. An extended Heisenberg model including a sizable 
four-site cyclic exchange term is shown to be consistent with the experimental data. 

Pacs Numbers: 78.30.Hv. 75.40.Gb. 75.40.Mg. 75.50.Ee 



The starting point of many theories describing high- 
temperature superconducting cuprates is the undoped 
parent compound. A consistent description of this phase 
is of great importance since the usual approach is to "ex- 
tend" this model to describe the doped phase like the t-J 
model 01 or the spin fermion model Q| . Also the under- 
standing of this phase, namely the physical realization of 
a two-dimensional (2d) spin-1/2 quantum antiferromag- 
net, is a fundamental problem in itself. 

It is usually assumed that the 2d Heisenberg model 
(HM) 
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describes the physics of the stoichiometric materials ||] 
with Jij — J for the first nearest neighbors and zero oth- 
erwise. 

In this work we show that infrared (IR) optical absorp- 
tion spectra due to phonon assisted multimagnon excita- 
tions is very sensitive to the magnetic Hamiltonian. We 
find that the usual Heisenberg description is incompatible 
with IR experimental data when more than two magnons 
are involved. An extended Heisenberg model including 
further neighbor interactions in Jij and a four-spin cyclic- 
exchange (4SCE) term [p|-[l0[ is shown to explain the IR 
experiments. The same model has been shown to be con- 
sistent with other spectroscopic data | |ic| ]. 

The 4SCE was introduced by Takahashi |^ and by 
Roger and Delrieu on the present context, using a 4th or- 
der perturbative analysis |^ . It was also supported by ex- 
act diagonalization (ED) studies of a multiband-Hubbard 
model 1^ describing the Cu-O planes. The main effect 
of this term in the Hamiltonian is to permute cyclically 
four spins on a plaquette. 

Though small exchange interactions, going beyond the 
first neighbors are expected pdl], the 4SCE term is not 



generally accepted in the literature Up to now there 
has not been any clear experimental way to rule out or 
confirm the presence of this term. In fact the experimen- 
tally established |^ single-magnon spin-wave branch is 
rather insensitive, at low energies, to the presence of the 
4SCE term 

Any noticeable effect of the 4SCE term is expected to 
occur at high energies which are accessible in optical ex- 
periments. One example is magnetic Raman (MR) light- 
scattering p^ . The line shape has an asymmetric peak 
close to 3 J due to two-magnon excitations and a shoul- 
der at higher energy which is believed to arise from a 
four-magnon process 0. The width of the two -magnon 
peak and the four magnon shoulder are anomalous in 
the sense that they do not agree with a conventional in- 
teracting spin-wave theory description of the line shape 
|p|- pl] , p^ 16 1 . Theoretical studies have attributed both 
anomalies to the presence of other terms in the Hamilto- 
nian §-|l|,|l|,|l6| including the 4SCE term §-0. Al- 
though the latter assignment is encouraging we argue 
that the analysis of the MR line shape is not conclusive 
(see Ref. @)- 

Another experiment, which probes the multimagnon 
response, is phonon-assisted multimagnon light absorp- 
tion (PAMLA) ||r^-|2^. In this experiment an absorbed 
photon simultaneously creates a phonon and a multi- 
magnon excitation. The absorption mechanism is well 
understood [^,^ allowing to make theoretical predic- 
tions on the nickelates p^(b ) (2d, spin 1) which were 
successfully corroborated [ |19pO| ]. Besides, the line shape 
was computed for spin-1/2 Id Cu-O systems |23j with 
great accuracy p2|(d). The experimental line shape is 
reproducible even among different materials (see Fig.^). 
All this puts the interpretation of the IR line shape on 
a firmer basis than the MR line shape |0] making it an 
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ideal candidate to test models of the spin dynamics. 

For the 2d cuprates the line shape has been measured 
in several materials p8|-^ . In Fig. [l] we show the line 
shape for three different cuprate compounds in a dimen- 
sionless scale. All the data collapse to a unique curve, 
implying that this experiment depends only on the com- 
mon CuO 2d layers. The high energy upturn is due to 
the charge transfer band and the low energy upturn is 
due to the phonons. 

The line shape has, like in MR scattering, a structure 
close to 3J (measured from the phonon energy) due to 
the two-magnon process. In addition strong side-bands 
appear at high energies. 

A recent study has suggested that, though the main 
peak is of magnetic origin, the side bands may be ex- 
plained by the presence of a d-d exciton ||2^. In this 
case, since a d-d exciton will depend on details outside 
the Cu-0 planes one would expect the position and in- 
tensity of the side bands to be unrelated to the position 
and intensity of the main peak when different materials 
are examined. Instead the scaling shown in Fig. ^ shows 
that this is not the case ruling out the exciton expla- 
nation. We anticipate that the main peak and the side 
bands can be understood within the PAMLA mechanism 
if the appropriate magnetic Hamiltonian is used. 
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FIG. 1. I{lj~ujo) obtained from the experimental absorp- 
tion using Eq. . Following Refs. [ p^ , po| a linear background 
has been subtracted form the raw absorption data. We show 
the data for La2Cu04 (l|, Sr2Cu02Cl2 and YBa2Cu306 
[ pl[ in a dimensionless scale using the following reference en- 
ergies (ojo, J) = (0.08,0.12), (0.06,0.11) and (0.06,0.1) eV re- 
spectively. The energy scale J and the intensity scale were 
adjusted to match the first peak in ED. We also show the 
HM theoretical line shape in ED (thin line) for a 32 site clus- 
ter and a Lorentzian broadening of 0.16 J and in interacting 
spin- wave theory j22j (dashed line) . The intensity scale of the 
latter was adjusted to match the ED intensity. 

The PAMLA absorption coefficient is given by, 



where ao is a material dependent constant defined in 
Ref. LjQ is the frequency of the stretching mode 

phonon. /(w) is the weighted sum of a two-spin spec- 
tral function on the Brillouin zone [pll 



" TV ^ sin(fc/2)2(sin(fc/2)2 + smiqy/2y 
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we introduced the Fourier transform = 



and 

_^ e^q-Rj Sr. . 

The two-magnon peak has been accurately fitted with 
a two-magnon interacting spin-wave theory computation 
|^(a),(b) in the HM which we report on Fig. |l|. 

The physics of the two-magnon peak and the side 
bands can be understood with a simple argument [^,^ . 
If one approximates the ground state by the classical Neel 
state the effect of the operator is to flip two spins in 
nearest neighbor sites. The energy of this excitation is 
3 J which is close to the energy of the peak observed. In 
the HM this state is not an eigenstate and will mix with 
states with four, six, etc. spin flips. Consequently the 
spectral function will show side bands at the energy of 
these excitations which in the Ising limit are 4 J, 5 J, etc. 

Since the side bands involve more than two magnons 
they cannot be described in the two-magnon-interacting 
spin- wave theory computation of Ref. |^(a),(b). It is 
therefore important to avoid the spin-wave approxima- 
tion and establish with an unbiased technique whether 
or not the side bands can be described by the HM. To 
this purpose we have computed the spectrum using ED 
on finite clusters. 

In Figs. we show the exact spectrum in different 
size clusters. Although the Heisenberg line shape has 
some structure at the energy of the side bands the rela- 
tive intensity is much smaller than the experimental one. 
Notice the similarity with the two-magnon spin- wave the- 
ory line shape (dashed line in Fig. |l|). 

To make a more quantitative comparison we have com- 
puted the cumulants |ll|, 
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with It = J I{u!)duj and p = ^ ujI{uj)dio j It- 

M2 and measure the width and asymmetry of I{ijj) 
respectively {Mi = 0). They are are obviously very sen- 
sitive to the presence of the side bands so it is natural to 
use them to characterize the line shape. 

In Fig. ^ we show M2 and and the ratio of the av- 
erage energy p to the same quantity in MR (pRaman) for 
different system size and for the experiment. 
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An extrapolation to an infinite system confirms what 
Figs. dJI suggest: the HM alone can not correctly de- 
scribe this experiment. It is difficult to ascribe this to a 
failure of the PAMLA mechanism itself since, as we men- 
tion above, the mechanism has been successfully tested in 
an isostructural system with spin one, namely La2Ni04 
[^(b) and a system with the same spin (1/2) but 
lower dimensionality (d) 1 23 1 . We therefore analyze a 
more realistic Hamiltonian to describe the spin dynamics. 




limited size cluster that the agreement with the experi- 
mental data is improved. As a consequence of this trans- 
fer of spectral weight one sees an increase of moments 
(Fig. I). 
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FIG. 2. I{uj) in ED for a 26 site cluster and with the same 
parameter sets as in Fig. ^. We also show the experimental 
line shape in Sr2Cu02Cl2 ^] (thick line). The intensities 
were scaled to coincide at the maximum of the HM spectrum. 
In order to have the maximum of all the spectra at the same 
UJ we have used different values of J for the three different 
models: J = O.leV (HM) J = 0.12 eV (4SCE with K/J=0.3) 
J = 0.15 eV (EHM) 

The Hamiltonian with the 4SCE term reads H = 
Huci + H^s with 

H4S=K iSi.SjKSk.Si) + iSi.Si){Sj.Sk) 

<i,j,k,l> 

-{Si.Sk){Sj.Si) (5) 

where < i, j,k,l > stands for the sum over groups of four 
spins on a plaquette. This term can be shown to produce 
the cyclic permutation of the four spins on the plaquette 
plus ordinary two spins exchanges of all the pair of spins 
of the plaquette including the ones on the diagonals (see 
Ref. [PH25I]). The parameter K/J has been estimated us- 
ing an ED mapping from a multiband Hubbard model j^] 
to be around 0.3. 

In Fig. H we show the exact line shape in the HM and 
in the model with the 4SCE term in a 26 site cluster. 

We notice a strong sensitivity of the spectra to the 
4SCE term. The main effect is to transfer weight from 
the first peak to the side bands. We already see in this 



FIG. 3. a) and b) Mj/p as a function of the inverse system 
size, c) p/pRaman [M • We show ED data for the HM (trian- 
gles), for the HM with 4SCE and K/J = 0.3 (circles) and in an 
EHM with K/J = 0.34, J2/J = 0.04, J3/J = 0.03 (squares). 
The dashed line is a linear 1/A^ extrapolation for A'^ > 16 in 
the HM. The dotted lines are guides to the eye traced paral- 
lel to the Heisenberg extrapolation. The experimental values 
are computed from the Sr2Cu02Cl2 data p^p9[ |. The error 
bar reflects indeterminacies in background subtraction and, 
for pRaman, the dependence on excitation energy |29| . 

A fourth order perturbative analysis of the Hubbard 
model jl^ generates further nearest neighbor exchange 
interactions, beyond the 4SCE term already considered. 
The additional terms are spin exchange interactions for 
next-nearest neighbor sites (Jjj — J2) and for next-next- 
nearest neighbor sites ( = J3) in Eq. (|]). In order to be 
systematic we consider therefore an extended Heisenberg 
model with all magnetic interactions arising at fourth or- 
der. 

To estimate the parameters we followed Ref. and 
apply the perturbative expressions to an extended Hub- 
bard model including second {t2nd/tist = 0.15) and third 
(tzrd/tist = —0.12) neighbor hoppings t. The Hubbard 
U is taken as Sti^j. This approach has been shown to be 
consistent with MR, neutron-derived spin wave velocity 
as well as angular-resolved-photoemission spectroscopy 
data ||l^. From now on we refer to the resulting mag- 
netic Hamiltonian as extended Heisenberg model (EHM). 

In Fig. ^ we show the line shape and in Fig. ^ we show 
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the size dependence of the moments. We see that the 
effect of these terms also improves the agreement. Again 
we notice a strong sensitivity of the spectra to the un- 
derlying magnetic Hamiltonian. 

A fine tuning of the Hamiltonian parameters, which ac- 
curately reproduce the experiments, is not possible with 
the limited sizes available, due to the difficulty of a pre- 
cise extrapolation of the cumulants to an infinite system. 
Clearly a sizable value of K/ J is needed and probably 
also a non negligible J2/J and J3/J. Our best estimate 
is the latter parameter set considered. In fact a rough 
extrapolation from Fig. assuming a similar scaling as 
for the HM, gives a result quite close to the experimental 
data. We mention that by setting K = Q m the EHM we 
where not able to obtain an acceptable fit to the exper- 
imental data with reasonable values of J2/J and J3/J. 
Even for the EHM Fig. || does not show perfect agree- 
ment between theory and experiment but the moment 
analysis shows that this can be ascribed to a finite size 
effect. 

The energy scale can be fixed by matching the first mo- 
ment with the experimental first moment. In the EHM 
we find J ^ Q.l'deV for Sr2Cu02Cl2 in good agreement 
with the value found in Ref. using other spectro- 
scopic data. As an alternative procedure one can adjust 
the position of the first peak as done in Fig. ||. This 
gives a somewhat smaller value probably due to finite 
size effects. These values of J are not in contradiction 
with the smaller value of J usually quoted in cuprates 
(J ^ O.I ~ 0.13 eV): A spin-wave theory computation 
shows that, at low energies, the effect of the extended 
terms in the Hamiltonian is to renormalize the effective 
J to lower values |^]. 

We also computed the staggered magnetic moment in 
the EHM. We get a staggered moment roughly 7% larger 
than in the HM (see also Refs. j|,|2^). This may be im- 
portant in view of the disagreement found for this quan- 
tity between theory and experiment |26y27f| . 

In conclusion we have presented a computation of the 
IR absorption spectra due to magnetic excitations in un- 
doped cuprates. We have shown that these experiments 
are very sensitive to the underlying magnetic Hamilto- 
nian. We find that the usual model used to describe the 
spin dynamics in cuprates, namely a HM with nearest- 
neighbor exchange, can not explain the experimental 
data. Instead an EHM with further neighbor interactions 
and with a 4SCE term is in good agreement with the 
data. To the best of our knowledge this provides the first 
quantitative explanation of the puzzling side bands in the 
spectrum. In addition, the same model has been found 
to agree with other spectroscopic data ||l^ . Further the- 
oretical work is needed to explore the consequences of the 
4SCE term in the doped phase. 

We have been strongly influenced by ideas and views 
from Eduardo Gagliano who unfortunately passed away 
during the course of this work. We are in debt with the 
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